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Relatively general techniques for computing mean first-passage time (MFPT) of random walks on 
networks with a specific property are very useful, since a universal method for calculating MFPT 
on general graphs is not available because of their complexity and diversity. In this paper, we 
present techniques for explicitly determining the partial mean first-passage time (PMFPT), i.e., the 
average of MFPTs to a given target averaged over all possible starting positions, and the entire 
mean first-passage time (EMFPT), which is the average of MFPTs over all pairs of nodes on regular 
treelike fractals. We describe the processes with a family of regular fractals with treelike structure. 
The proposed fractals include the T fractal and the Peano basin fractal as their special cases. We 
provide a formula for MFPT between two directly connected nodes in general trees on the basis of 
which we derive an exact expression for PMFPT to the central node in the fractals. Moreover, we 
give a technique for calculating EMFPT, which is based on the relationship between characteristic 
polynomials of the fractals at different generations and avoids the computation of eigenvalues of 
the characteristic polynomials. Making use of the proposed methods, we obtain analytically the 
closed-form solutions to PMFPT and EMFPT on the fractals and show how they scale with the 
number of nodes. In addition, to exhibit the generality of our methods, we also apply them to the 
Vicsek fractals and the iterative scale-free fractal tree and recover the results previously obtained. 

PACS numbers: 05.40.Fb, 61.43.Hv, 89.75.Hc, 05.60.Cd 



I. INTRODUCTION 

As an integral ingredient of nonlinear science, frac- 
tals have attracted an increasing attention in physics and 
other scientific fields [l], H| not only because of the strik- 
ing beauty intrinsic in their structure Q but also due to 
the significant impact of the idea of fractals on a large va- 
riety of scientific disciplines, such as astrophysics, plasma 
physics, optics, economy, ecology, and so on |4{. Among 
various fractal categories, regular fractals constitute an 
important family of fractals. Frequently cited examples 
include the Cantor set [f|, the Koch curve @, the Sier- 
pinski gasket 0] , the Vicsek fractals [42j , the T fractal [t| , 
etc. These structures have been a focus of research ob- 
ject [l(| • One of the main justifications for studying regu- 
lar fractals is that many problems can be exactly solvable 
on regular fractals [111 ], the explicit solutions for which 
provide useful insight different from that of approximate 
solutions for random fractals. 

One of the ultimate goals of research efforts on frac- 
tals is to unveil how their underlying geometrical and 
structural features affect dynamical processes and criti- 
cal phenomenon on them [12j . Among a variety of dy- 
namics, due to their simplicity and wide range of appli- 
cations (l3| . random walks play a central role in many 
branches of sciences and engineering and have gained 
a considerable attention within the scientific commu- 
nity [li], EH- A key quantity related to random walks 
is mean first-passage time (MFPT) [l6|, [TtJ that is the 



expected time for the walker starting off from a source 
node to first reach a given target node. In addition to its 
role as a basic quantitative measure of the transporta- 
tion efficiency, MFPT also encodes useful information of 
other quantities concerned with random walks [l6j . 

Concerted efforts have been devoted to study MFPT in 
fractals in order to obtain the scaling of MFPT with sys- 
tem size [l8| . For instance, MFPT associated to a given 
starting node has been discussed by many authors [l9l . l20j 
since it depends on the starting position and is an impor- 
tant issue in its own right. On the other hand, recent pa- 
pers have addressed the explicit determination of partial 
mean first-passage time (PMFPT) on some regular frac- 
tals, defined as the average of MFPTs to a selected tar- 
get node averaged over allpossible starting nodes, such 
as the Sierpinski gasket [H H, the T fractal HE HI, 
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the iterative scale-free treelike fractal [25[, as well as the 
hierarchical lattice fractals (2(|[27]]. Besides, entire mean 
first-passage time (EMFPT), i.e., the average of MFPTs 
over all pairs of nodes, has been computed for the T frac- 
tal [HI, the Vicsek fractals (2^, and some other regular 
fractals [Hj]. Thus far, most techniques used to calcu- 
late PMFPT or EMFPT are only applicable to specific 
structure; universal (even relativel y g eneral) methods for 
computing MFPTs are much less jl7l |30| . 

In this paper, we develop techniques for the explicit 
determination of PMFPT and EMFPT on regular tree- 
like fractals, which are also expected to be valid to other 
deterministic media. For PMFPT, we provide a univer- 
sal formula for MFPT between two adjacent nodes in 
generic trees based on which one can determine PMFPT 
between a given node and all other nodes. For EMFPT, 
by considering only the relationship between characteris- 
tic polynomials at different iterations, we can find exactly 
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FIG. 1: (Color online) Iterative construction method of the 
treelike regular fractals. The next generation is obtained by 
performing the operation shown to the right of the arrow. 



the analytical expression for EMFPT, which needs only 
local information of characteristic polynomials and can 
avoid the laborious explicit computation of the eigenval- 
ues. 

To illustrate the calculation processes, we put forward 
a family of treelike regular fractals, which include the fa- 
mous T fractal and the Peano basin fractal [3lJ as their 
particular cases. Using the presented methods, we first 
study the random-walk problem with a trap fixed at the 
central node of the fractals and obtain the PMFPT to 
the trap; then, we address the case that the trap is uni- 
formly selected among all nodes and derive the rigorous 
solution to EMFPT between all node pairs. From the 
obtained analytical results, we give the scalings for both 
PMFPT and EMFPT and show that they both increase 
as a power-law function of the number of nodes, with the 
exponent larger than 1 but smaller than 2. 



II. CONSTRUCTIONS AND PROPERTIES OF 
THE TREELIKE FRACTALS 

Let us first introduce the model for the treelike frac- 
tals that are built in an iterative way and controlled by 
a positive-integer parameter m (i.e., m > 1). We de- 
note by T g (g > 0) the treelike fractals after g iterations. 
Notice that the number of iterations also represents the 
generation of the treelike fractals. Initially (g = 0), To is 
an edge connecting two nodes. For g > 1, T g is obtained 
from Tg-i by performing the following operations on ev- 
ery edge in T g _\ as shown in Fig. [T] replace the edge by 
a path of two links long, with the two end points of the 
path being the same end points of the original edge, then 
attach m new nodes to the middle node of the path. Fig- 
ure [2] illustrates the construction process of a particular 
fractal for the case of m = 2, showing the first several 
generations. Notice that there are several limiting cases 
of our model. When m = 1, it is reduced to the T frac- 
tal [9j. When m — 2, it turns out to be the Peano basin 
fractal 

According to the construction algorithm, at each step 
the number of edges in the treelike fractals increases by 
a factor of m + 2. Thus, we can easily know that the 
total number of edges of T g is E g — (m + 2) 9 and that 
the total number of nodes (often named network order) 
in T g is N g = E g + 1 = (m + 2) 9 + 1. One can partition 
all the N g nodes into two classes: nodes having degree 
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FIG. 2: Illustration for the growing process for a special frac- 
tal corresponding to the case of m = 2. A circle expresses an 
external node, while a square denotes an internal node. 

m + 2 are called internal nodes and nodes with degree of 
1 are named external nodes. Since at each generation, 
a preexisting edge creates m + 1 nodes: one of which 
has a degree of m + 2, the rest m nodes have a degree 
of 1. Thus, at generation gi (gi > 1), the numbers of 
internal and external nodes generated at this generation 
are (m + 2) 9i ~ 1 and m(m + 2) 9i ~ 1 , respectively. 

On the other hand, it is obvious that after each it- 
eration, the diameter (namely, the maximum of short- 
est distances between all pairs of nodes) of the fractals 
doubles. Thus, we have that the fractal dimension of 
the treelike fractals is dj = ln(m + 2)/ In 2. In addi- 
tion, as we will show below, for any two nodes i and j 
at current generation, the MFPT from node i to node 
j increases by a factor of 2(m + 2) at next generation. 
Hence, the random-walk dimension of the fractals is d w — 
ln[2(m + 2)]/ In 2 = 1 + df, and their spectral dimension 
is d= 2d f /d w = 21n(m+2)/ln[2(m+2)] = 2d f /(l + d f ). 

It is worth mentioning that the treelike fractals can 
also be alternatively constructed using another method. 
If we define the central node (see Fig. [5] for an example) 
as the innermost node and the nodes farthest from the 
central node as the outermost nodes, then the second 
generating algorithm for the fractals can be described as 
follows: Given the generation g, T g+ i can be obtained 
by joining m + 2 replicas of T g (see Fig. [3]) . That is to 
say, to obtain T g+ i one can merge together the discrete 
outermost nodes of the m + 2 copies of T g . The rn + 2 
outermost nodes belonging to separate copies merge into 
a single new node, which is then the innermost node in 

Tg + l- 

After introducing the construction algorithms and the 
properties of the treelike fractals, in what follows we will 
study the random-walk dynamics on the fractal family 
Tg. First, we will investigate random walks with a single 
immobile trap located on the innermost node. Then, 
we will continue to address random walks with the trap 
selected randomly (i.e., distributed uniformly) among all 
nodes. 
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FIG. 3: (Color online) An alternative construction method of 
the treelike fractals. The fractals after g + 1 iterations, T a +i, 
can be obtained by merging m + 2 copies of T g , sequentially 



denoted by T 9 (1) , T< 



.(2) 



T, 



(m + l) 



and T, 



(m + 2) 



where X is 



the innermost node, while K is an outermost node. 



III. RANDOM WALKS WITH A TRAP FIXED 
ON THE CENTRAL NODE 



In this section we focus on a particular case of random 
walks — trapping problem — on the family of treelike frac- 
tals T g with the trap or perfect absorber placed on the 
central node of the regular fractals. In the process of the 
random walks, at each time step, the walker, starting 
from its current location, jumps to any of its neighbor 
nodes with identical probability. For the convenience of 
description, we label by 1 the central node of T g , while all 
other nodes are labeled consecutively as 2, 3, . . ., N g — 1, 
and N g . Let Fi(g) be the MFPT, also called trapping 
time, for a walker, staring from node i to first arrive at 
the trap. What we are concerned is the PMFPT pre- 
sented as (F) g , which is defined as the average of Fi{g) 
over all starting nodes distributed uniformly over nodes 
in the proposed fractals other than the trap. PMFPT 
is also often called mean trapping time or mean time to 
absorption 23[. By definition, (F) g is given by 



(F) 3 = 



1 



(i) 



i=2 



The main purpose of this section is to determine explic- 
itly (F) g and show how (F) g scales with network order. 



A. Mean first-passage time between two adjacent 
nodes in general trees 

Here, we provide a universal formula for MFPT be- 
tween two nodes directly connected by an edge in a 
treelike network, which is very helpful for the following 
derivation of (F) g . Let us consider a general connected 
tree. We denote by e = (u, v) an edge in the tree connect- 
ing two nodes u and v. Obviously, if we remove the edge 
(u, v), the tree will be divided into two small subtrees: 



one subtree includes node u and the other contains node 
v. We use C u<v to denote the number of nodes in the 
subtree containing node u. Actually, C u<v is the num- 
ber of nodes in the original tree lying closer to u than 
to v, including u itself. Furthermore, let F uv denote the 
MFPT from u to v. Then, F uv can be written in terms 
of Cu<v as 



— 2C- U < V 1. 



(2) 



Equation ^ can be easily proved inductively. If it is 
a leaf node (namely, node with degree 1), it is obvious 
that C u<v = 1 and F uv = 1, and thus Eq. ([2]) holds. We 
proceed to consider the case that u is a non-leaf node. For 
this case, to prove the validity of Eq. ©, we consider the 
tree as a rooted one with node v being its root, then v is 
the father of u and C u<v is in fact the number of nodes in 
the subtree whose root is u. Suppose that for a non-leaf 
node u, Eq. ^ is true for all its children, the set of which 
is denoted by f2 u . In other words, for an arbitrary node 
x in Q u , the following relation holds: 



(3) 



Then, the MFPT F MV from node u to its father node v 
can be calculated by 



F u 



1 1 



£ (1 + Fxv) , 



(4) 



where 6 is the degree of node u. 

The first term on the right-hand side (rhs) of Eq. Q 
explains the case that the walker, starting off from node 
u, jumps directly to node v in one single step with proba- 
bility | . The second term accounts for another case that 
the walker first reaches a child node x of node u in one 
step, and then it takes more F xv steps to first hit the 
target node v. Equation (fj| can be rewritten as 



(5) 



which leads to 



J- II, 



l+£(! 

2C ll <„ — 1, 



F 



y ' (2C K<ll — 1) 

i6f!„ 



(6) 



where the assumption, viz. Eq. ([3]), was used. Thus, we 
have proved Eq. ([2]), which is a basic characteristic for 
random walks on a tree and is useful for the following 
computation of the key quantity (F) g . 



Evolution law for mean first-passage time in the 
fractals 



Let Fij (g) be the MFPT of random walks on T g , star- 
ing from node i, to first reach node j. Let [u, v) be 
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FIG. 4: Illustration for the evolution of the fractals. 

an edge connecting two nodes u and v in T g . Next, we 
will give a relation between F uv (g + 1) and F uv (g). To 
this end, we look upon the treelike regular fractals T g as 
rooted trees with node v as the root, and thus v is the 
father of u. Note that in the evolution of the fractals, 
we suppose that node v is always the root. In addition, 
for T g we use C u (g) and E u (g) to represent the numbers 
of nodes and edges in the subtree, whose root is node u. 
Clearly, 

E u (g) = C u (g) - 1 . (7) 

On the other hand, by construction of the fractals, each 
edge in T g will be replaced by to+2 new edges at iteration 
g + 1 . Then, we can easily have 

E u (g + l) = (m + 2)E u (g), (8) 

from which we can obtain the recursion relation between 
Cu(g + 1) and C u (g). By definition, 

C u {g + l) = E u (g + l) + l = (m + 2)E u (g) + l 

= (m + 2)C u {g)-m-l. (9) 

Now we begin to derive the relation governing F uv (g + 
1) and F uv (g). According to the general result given in 
Eq. ([2]), we have 

F vv {g) = 2C u {g)-l. (10) 

Notice that at generation g+1, v is no longer the father of 
u but still an ancestor of u. According to the construction 
algorithm of the fractals, the edge in T g connecting nodes 
u and v will generate to + 1 new nodes at generation 
g + 1 (see Fig. [5]): one node x has a degree to + 2, while 
the other to nodes denoted by Xi (i = 1, 2, . . . , to) have 
the same degree 1. After this iteration, the father of u 
becomes x and a: is a child of v. Then, for a random 
walker in T g+ i, if it wants to transfer from u to v, it 
must pass by node x. Therefore, 

F uv {g + 1) = F ux {g + 1) + F xv (g + 1) 

= 2C u (g + l)+2C x (g + l)-2. (11) 



Then, the determination of F uv (g + 1) is reduced to find- 
ing C u {g + 1) and C x (g + 1). 

Since C u (g + 1) has been determined above, we only 
need to find C x (g+ 1). According to the structure of the 
fractals, it is obvious that 

C x {g + 1) = C u (g + 1) + to + 1 = (to + 2)C u {g) , (12) 

where the term to + 1 between the two equal marks de- 
scribes the number of nodes x and xi (i — 1,2,..., to). 
Plugging Eqs. © and (H]) into Eq. dTTJ) , we obtain 

F uv (g + 1) = 4(TO + 2)C u ( ff )-2TO-4 

= 2(m + 2)F uv (g), (13) 

in which Eq. (jTDJ) has been made use of. 

Equation (fT3")) tells us that for any two adjacent nodes 
u and v in the fractals at a given generation, the MFPT 
from one of them to the other will increase by a factor 
of 2(to + 2) at the next generation. Since the consid- 
ered fractals have a treelike structure, for any pair of two 
nodes i and j, adjacent or not, the MFPT between them 
obeys the relation: 

F ij (g + l) = 2(m + 2)F ij (g), (14) 

which is a basic property of random walks in the regu- 
lar treelike fractals, dominating the evolution of MFPT 
between any couple of nodes. The relation described by 
Eq. (|14[) is very important using which we will derive the 
rigorous formula for {F) g . 

C. Exact solution to partial mean first-passage 
time in T g 

Having obtained the evolution rule of MFPT for ran- 
dom walks in the fractals, we now determine the mean 
time to absorption averaged over all non-trap nodes in 
T g . To attain this goal, we represent the set of nodes in 
T g as A g , and denote the set of those nodes created at 
generation g by A g . Obviously, A g = A g _i + A s . For the 
convenience of description, we define the following two 
quantities for n < g: 

F^(g)=J2 F *(9), (15) 

ieA„ 

and 

F^(g)= (16) 

ieA„ 

Thus, we have 

Fl°\g) = 2(m + 2)Ff\{g) + F^( 9 ) , (17) 

where we have used Eq. ((14)) . Thus, to obtain Fg 0t (g), 
we should first determine the quantity Fg Qt (g). 
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To find Fg 0t (g), we separate set A g of the nodes created 
at generation g into two subsets A g nt and A° xt , such that 
A g = Ag Xt U A g nt , where A g nt is the set of internal nodes 
and A g xt is the set of external nodes. As shown in Sec. 
III1 the cardinalities (the cardinality of a set is the number 
of nodes in the set) of the two subsets, denoted by |A® xt | 
and |A g nt |, are |A g nt | = [m + 2) 9 - 1 and |A g xt | = m(m + 
2) 9_1 , respectively. Clearly the two variables satisfy the 
following relation: 



\K 



m|A 



int 



Then, the quantity F g ot (g) can be rewritten as 

^ tot (s)= E w+ E 



(19) 



which shows that to determine F* ot (g), one may alter- 
natively find the two quantities on the rhs of the equal 
mark. We begin from determining the first sum term 

T,i E A^ F i(9)- 

Since at a given generation g, [gi > 2), any internal 
node has a degree of m + 2 and two internal nodes in 
A™* are not directly connected to each other, we have 
that any of the E gi edges in network T ai must have one 
(and only one) end point belonging to A^'. In addition, 
by construction any internal node in is linked to 

both ends of the edge preexisting at generation gi that 
generated the internal node, which means that any inter- 
nal node in Aj^^i m ust have one (and only one) neighbor 
belonging to A™'. Furthermore, we can know that at gen- 
eration gi + X, the m + 2 neighbors of any internal node 
in A™* are all internal nodes in A™+ ls and for any two 
internal nodes in A™*, their neighbors are completely dif- 
ferent. From this information, we can obtain a recursion 
relation for J2ieA™< F i(d)- 

Consider the random walks on the fractals T g . For an 
internal node i 1 g a ^ 1 in A^ 1 and all its m + 2 neighbors 
1, 2, . . . , m + 2) belonging to A ff nt , their trapping 
times obey the following relation: 



F^(g) = X + 



1 



m+2 



(20) 



Summing Eq. (f20|) over all internal nodes in A^j, we 
obtain 

E m = \x£x\ + ^T2 S (21) 



i£A» 



which may be rewritten as 
E Fi(g) = (m + 2) E F l {g)-{m + 2)\KfU\ 

= 2 (m + 2) 2 E F i {g-X)-{m + 2) g - 1 . 



i£A" 



(22) 



Usin S EieAjf* ^( 2 ) = ( 2 ™ + 3)(m + 2), Eq. (22) is re- 
solved by induction, 

E F *(») = ( ^ + ?1 2 [ 2 + m + 29 (™ + + 2) 9 ] 



ieA" 



3 + 2m 



(23) 

After obtaining the first sum EieA lnt F i(g) in Eq. (fT9|) . 
we proceed to determine the second sum EiGA oxt Fi(g). 
A random walker staring from a node i e xt belonging to 
(18) A° xt must be on its only neighbor i- mt — an internal node 



Fi^(g) = F iint (g) + ^ 



(24) 



Summing Eq. ([24]) over all external nodes in A g gives 
rise to 



E F i (g)=m E W + 



CXt I 
5 I 



(25) 



ieA= 



iGA" 



where we have considered Eq. (|18|) . Substituting Eq. (|23|) 
for EieAj- and |A° xt | = m(m + 29- 1 ) into Eq. flU), 
we have 

Fi(g) = m( To! )9 (2 + 7n + 2 9 (m+l)(2 + m) g ) 



3 + 2m 
+m (m + 2) 9-1 



(26) 



With the results obtained for EieA int F i(d) an< ^ 

EiGA<= xt Fi(g), we can determine F g ot (g). Substituting 

the two sums given in Eqs. (l23l) and (j26|) for F g ot (g) into 
Eq. dT7J) yields to 

F^ig) = 2(m + 2)F 9 t _ t 1 (.g) + m(m + 2f- 1 + 
(m+ l)(2 + m) 9 ~ 2 



3 + 2m 



[2 + 771 + 2 s (m + 1)(2 + m) 9 ] . 

(27) 



Considering the initial condition F' ot (l) = m + 2, 
Eq. (|27|) is resolved by induction to yield 



F^ig) 



(2 + m) 9 



■[2 9+1 (m+l)(2 + m) 9 



9 vay 2(2m 2 + 7m + 6) 

+ (2m 2 + 5m + 3)2 9 - 4m 2 - 8m - 2] (28) 

Plugging the last expression into Eq. (QJ, we arrive at the 
explicit formula for the PMFPT on the studied treelike 
regular fractals, 



1 N ' 1 



9 i=2 
1 



[2 9+1 (m + X)(2 + m) B 



2(2m 2 + 7m + 6) 
+ (2m 2 + 5m + 3)2 9 - 4m 2 - 8m - 2] . (29) 
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When m = 1, Eg. (|29j) recovers the result previously ob- 
tained in Refs. [23[ and |23| . confirming that the solution 
given by Eq. (|29|) is right. 

We continue to show how to represent PMFPT as a 
function of network order N g , with the aim to obtain the 
scaling between these two quantities. Recalling N g — 
(m+2)s+l, we have g = log m+2 (iV a -l). Hence, Eq. p?]) 
can be recast as 

( F ) 9 = 0,0^2 hz. , g , [2(m + l)(JV g 



average: the first one is over all the starting nodes to a 
given trap node and the second one is the average of the 
first one over all trap nodes. By definition, (H) g is given 
by 



N„ 



(H) 



N g (N g - 1) ^ 



(32) 



- 7m + 6) 
t-5m + 3)(iV s - 1) 
8to - 2] . 

For systems with large order, i.e., N g 



2(2m 2 - 
+ (2m 2 
—Am? - 



1 



(F) g ~ (iV g )l+ln2/ln(m+2) = (jy ) 2 / d 



i+in2/in(m+2) In general, one can apply the method for pseudoinverse 
of the Laplacian matrix 34] to calculate numerically but 
in2/in(m+2) exactly (H) g . However, this approach makes heavy de- 

mands on time and computational resources and is thus 
not tractable for large networks. To overcome these 
shortcomings, we will assort to the theory between elec- 
trical networks and random walks to derive {H) g , pro- 
(31) viding a closed-form expression for this related quantity. 



(30) 



oo, 



This confirms the lower bound for PMFPT provided in 
for general graphs (27J- Thus, for the whole family of the 
treelike fractals, in the large g limit, the PMFPT grows 
as a power-law function of the network order with the 
exponent, represented by 0(m) — 1 + In 2/ ln(m + 2), be- 
ing a decreasing function of m. When to grows from 1 
to infinite, the exponent 9{m) drops from 1 + In 2/ In 3 
and approaches to 1, indicating that the PMFPT grows 
supcrlinearly with network order. This also means that 
the efficiency of the trapping process depends on param- 
eter to: the larger the value of to, the more efficient the 
trapping process. 

Before closing this section, we stress that although we 
only focus on a special family of fractals, above compu- 
tational method and process are also applicable to the 
trapping problem on other self-similar trees. For exam- 
ple, we have used this method to compute the PMFPT 
to a target hub node in an iterative scale-free tree first 
introduced in [32[ and studied in detail in [33[ and to re- 
produce the result on PMFPT previously reported 25 1. 



IV. RANDOM WALKS WITH THE TRAP 
DISTRIBUTED UNIFORMLY ON THE 
FRACTALS 

In the Sec. Ill C, we have studied the PMFPT from a 
node to a trap fixed on the central node of treelike fractals 
averaged over all possible staring points, and showed that 
the PMFPT scales algebraically with network order. In 
this section, we will investigate random walks in T g with 
the trap uniformly distributed among all nodes of the 
fractals. As will be shown, the EMFPT (H) g , which is 
the average of MFPTs over all pairs of nodes, exhibits 
the same scaling as that of (F) g . 



A. Mean first-passage times and effective 
resistance 

The quantity EMFPT (H ) g concerned in this case that 
trap is randomly selected from all nodes involves a double 



To obtain the solution to 



H) g , we view T g as resistor 



networks [35j by considering all edges of T g to be unit 
resistors. Let Rij(g) be the effective resistance between 
two nodes i and j in the electrical networks obtained from 
T g . Then, according to the connection between MFPTs 
and effective resistance [30, S3 , we have 



F ij (g)+F ji (g) = 2E g R ij (g), 
using which Eq. (|32l) can be recast as 

N a No 



(33) 



E„ 



W-T)E E ^) 

u u ' i=i ]=i,j=£% 

, N g N g 

= wT, E *«<9), 



(34) 



t=i j=ij^i 



where the sum of effective resistance between all pairs 
of nodes is the so-called Kirchhoff index '38], which we 
denote b y -R tot(ff)- Using the previously obtained re- 
sults (HjliO], i?tot(s) can be rewritten as 

Rtot(g) = J2 E RiM=2N g J2l-7Z-y (35) 

i=lj=l,j^i i=2 l{9) 

where Xi(g) (i = 2, . . . , N g ) are all the nonzero eigenval- 
ues of the Laplacian matrix corresponding to T g . Note 



that Eq. 



only holds for a tree. Then, the EMFPT 



(H) g can be expressed as 



N„ 



(H) 



= ^E 



(36) 



We use L g to represent the Laplacian matrix for T g . 
The entries Lij(g) of L g are defined as follows: the off- 
diagonal element Lij(g) is —1 if the pair of nodes i and j 
are directly linked to each other, otherwise Lij(g) equals 
0, while the diagonal entry Lu(g) is equal to the de- 
gree of node i. After reducing the computation of (H) g 
to determining the sum of the reciprocal of all nonzero 
eigenvalues of L g , the next step is to evaluate this sum. 
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B. Using eigenvalues of the Laplacian matrix to 
determine entire mean first-passage time 

Let P g {\) express the characteristic polynomial of ma- 
trix L g , i.e., 

P s (A)=det(L s -AI s ), (37) 

where I g is an N g x N g identity matrix. As mentioned 
above, one of the main goals is to find the sum of the 
reciprocal of all nonzero eigenvalues of L g , namely, all 
nonzero roots of polynomial P g (X). 

In order to determine this sum, we denote Q g as a 
(N g — 1) x (N g — 1) sub-matrix of (L g — Xl g ), which is 
obtained by removing from (L g — AI g ) the row and col- 
umn corresponding to an outermost node, e.g., node K in 
Fig. G3 Also, we define R g as a sub-matrix of (L g — Xl g ) 
with order (N g — 2) x (N g — 2), obtained by removing 
from (Jj g — Xlg) two rows and columns corresponding to 
two arbitrary outermost nodes in two different copies, 
Tf\ (i = 1, 2, . . . , m + 1, m + 2), of T g _ 1 that are con- 
stituents of T g (see Fig. [3]) . On the other hand, we denote 
Q g (X) and R g (X) as the determinants of Q g and R g , re- 
spectively. Then, the three quantities P g (A), Q g (X), and 
Rg(X) satisfy the following recursion relations: 



Pg+lW 



Qg+lW 



m + 2 - 


A 


-e. 


-e, 






Qg 
O 


O 

Q 9 






O 


o 


m + 2 - 


A 




-e, 






Rg 

O 


O 

Q 9 






O 


o 



and 



Rg+l(X) 



m + 2 - A 
-e r T 
-eJ 



R g O 
O R 



y 



O O 



-e q 
O 
O 



o 
o 



Qg 

-e, 

o 
o 

Qg 



(38) 



(39) 



(40) 



In Eqs. (|3"8")) - (|4T))) . the superscript T of a vector represents 
transpose and e q (e r ) is a vector of order N g — 1 (N g — 2) 
with only the first entry being 1 and all other N g — 2 
(N g — 3) entries equaling 0, i.e., 



(1,0,0,0,. ..,0,0) 



(41) 



N„-2 



In Appendix [K\ using the elementary matrix opera- 
tions we show that P g+ i(A), Q g+ i(A), and R g+ i(X) evolve 
as: 

P g+1 (X) = (m+2)[Qg(X)r+ 1 P g (X) + (m+l)X[Q g (X)r+ 2 , 

(43) 

Qg+iW = [Q g (X)r +2 + (m + l)XRg(X)[Q g (X)] m+1 
+(m + l)R g (X)[Q g (X)] m P g (X) (44) 

and 

R g+1 (X) = 2Rg(X)[Q g (X)] m+1 + (m + l)A[ J R 9 (A)] 2 [Q g (A)] m 
+ m [R g (X)f[Q g (X)r~ 1 P g (X). (45) 

After obtaining the recursion relations for P g (A), 
Q g (X), and R g (X), next we will compute the sum of the 
reciprocal of the non-zero roots of P g (X). Since P g (X) 
has one and only one root equal to zero, say A \(g) = 0, 
to find this sum, we define a new polynomial P' g (X) as 



1 



PgW = TPgW 



Obviously, we have 

A'., 



Na-1 



(46) 



(47) 



where X[(g), X' 2 (g), ■ . ■ , A^ Vg _ 1 (.g) are the N g — 1 roots of 
polynomial P' g (X). Then, we reduce the problem to find- 
ing the sum on the rhs of Eq. (|4"T|) . 

Note that one can also express the polynomial P' g (X) 

in the form of P' g (X) = E^o' PgU)^ i h( * c Pgti) is the 
coefficient of term A- 7 with degree j] , such that 

N g -i N g -i 

£ P'gijW =P'g(Ng 1) ]J [A - X' L ( 9 )]. (48) 
3=0 i=l 



Comparing the coefficients of both sides of Eq. (148)) . we 
have 



Na-l 



P'g(l) 
1 K(9) ~ P'giOY 



51 Tu 



(49) 



Thus, we turn our aim into determining the two coeffi- 
cients p' g (0) and P g (l)- 

From Eqs. (|4"3")) - (l4"5|) we can easily get the recursive 
equations, 



P g+ iW = (m + 2)[Q g (A)] m+1 P g (A) + (m+l)[Q g (A)] 



m+2 

(50) 



Q g+1 (X) = [Qg(X)} m + 2 + (m + l)XRg(X)[Qg(X)Y 



]m+l 

+(m + l)XR g (X)[Q g (X)] m P g (X), (51) 



and 



and 



(1,0,0,0,. ..,0,0). 



A', 



(42) i? g+ i(A) = 2 J R g (A)[Q g (A)] m+i + (m+l)A[ J R g (Ar[Q g (A)] m 
+mA[i? g (A)] 2 [Q g (A)] m - 1 P g (A). (52) 
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Based on these relations we can find the values for p' g (0) 
andp s (l). 

We first evaluate p' g (0). For this purpose, let q g (0) 
and r g (0) be the constant terms of Q g (X) and R g (X), 
respectively. According to Eqs. ([5"U |) - (f52"1) . the quantities 
p' g (0), q g (0), and r g (0) can be represented recursively as 
follows: 

fWO) = (m + 2)[ <Z9 (0)] m +X(°) + (™ + !)fe(0)] m+2 , 

(53) 



m+2 



and 



%+i(0) - fe(0)] 



r 9+1 (0) = 2r g (0)< +1 (0). 



(54) 



(55) 



Using the initial conditions p'i(0) = — m — 3, gi(0) = f 
and ri(0) = 2, Eqs. (|53 ]) -([55 |) can be solved to yield 



p'(0) = -(m + 2)s -f 



9,(0) = 1 , 



(56) 



(57) 



and 



r g {0) = 2". (58) 

With these obtained results, we go on to compute 
p' g {X)- To this end, let q g (l) denote the coefficient of 
term A, i.e., term with degree 1, of polynomial Q g (X). 
Then, using Eqs. ([50 ]) - ([52 ]) . we have 



P ' g+1 (l) = (m + 2)[fc(0)P+y ff (l) 

+ ( m + 2)(m + l)[q g (0)} m q g (l)p' g (0) 

+ (m + 2)(m + l)[q g (0)} m+1 q g (l) (59) 



and 



q g+1 (l) = (m + 2)[q g (0)r +1 q g (l) 

+ (m + l)r 3 (0)[ (Z9 (0)r +1 

+(m + l)r 9 (0)[g g (0)]> a (0). (60) 

Substituting Eqs. ([56 > (|58" ]) into Eqs. flSSJ) and (JBDJl, we 
obtain the explicit expressions for p' g {l) and <z 9 (l) as 



9-1 



random walks in T g as 



(TO + 2) ~ [2 ff (m + 2)f+ 2 - 2(29 _ 1)(m + 2 )«rt-i 
2m + 3 

+ (2 9 -l)(m + 2) 9 + (m + 2) 2 -m-2l, (61) 



can arrive at the explicit formula of the EMFPT (H) g for 

2(m + 2)9~ 1 



N„-l 



^ Aj(ff) ((m + 2)9 + l)(2m + 3) 
[2 9 (m + 2) 9+2 - 2(2 5 - l)(m + 2) 9+l 
+ (2 9 - l)(m + 2) 9 + (m + 2) 2 - m - 2](.63) 



We have checked this rigorous result given by Eq. (|63|) 
against that obtained via numerical calculations by us- 
ing the method of pseudoinverse matrix [34| , which fully 
agree with each other. 

Similar to (F) g , we can rewrite EMFPT (H) g as a 
function of network order N g as 



(H) 



2{N g - 1) 



9 ' (2m 2 + 7m + 6)N g 

[(m 2 + 2m + l)(JV fl - i)i+in2/m (m +2) 

+ (2m + 3)(N g - 1) +m 2 + 3m + 2] . (64) 



Therefore, in the thermodynamic limit (N g — > oo), we 
have 

showing that the EMFPT (H) g grows as a power-law 
function of the network order N g . It should be mentioned 
that Eq. (|65p is in complete agreement with the general 
result for fractals given by Eq. (9) in [27] . 

Equations ([31]) and (|65[) indicate that (F) g and (i?) s in 
the treelike regular fractals show a similar behavior, both 
of which grow approximately as a power-law function 
of network order N g with the exponent 9{m) = 2/d = 
1 + In 2/ ln(m + 2) being a decreasing function of m and 
lying between 1 and 2, implying that they both increase 
superlinearly with network order. The sameness of scal- 
ings for (F) g and (H) g in T g is in comparison with the 
results obtained for the deterministic recursive tree [4lj |. 
where the PMFPT to a central node (a node with high- 
est degree) varies lineally with network order, while the 
EMFPT averaged over all node pairs grows with network 
order N as N\nN, a scaling larger than the linear one 
for PMFPT. 

Finally, we show the universality of our method for 
computing EMFPT in self-similar treelike fractals. Us- 
ing a process similar to that above, we determine the 
EMFPT in the Vicsek fractals, and recover the exact so- 
lution presented in one of our previous papers [!§]. For 
details, see Appendix IbI 



q g (l) = -(m + 2) 9 - 1 [1 - 2 s + 2 a (m + 2)] . (62) 

Thus far, we have obtained the intermediate quanti- 
ties, we now be in position to determine EMFPT (H) g . 
Inserting Eqs. ([5"6"]) and (|5T|) into Eq. (j3HJ), from which we 



V. CONCLUSIONS 

In this paper, we have presented techniques to deter- 
mine the explicit solutions to PMFPT and EMFPT in 
regular treelike fractals. First, we provided a formula 



9 



for the MFPT between two adjacent nodes in any con- 
nected tree, using which one can find PMFPT to a given 
target node. We then proposed a method for comput- 
ing EMFPT, i.e., average of MFPTs over all node pairs, 
which is based on the relationship of characteristic poly- 
nomials of the fractals at different iterations. It only 
needs partial knowledge of the polynomials, but avoids 
the laborious computation of eigenvalues that some other 
methods need. Thus, it is a simple and elegant method. 
We have used in detail our techniques to a family of tree- 
like regular fractals, which include some well-known frac- 
tals, e.g., the T fractal and the Peano basin fractal, as 
its peculiar cases. We also applied the methods to other 
treelike fractals to demonstrate their universality. 
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Appendix A: Derivation of recursion relations for 
relevant polynomials 

In this appendix, we give the detailed process for the 
derivation of recursive relations for polynomials P 9 (A), 
Q g (X), and R g (X). We first derive in detail the evolu- 
tion equation for P g (A), while the evolution relations for 
Q g {\), and R g (X) can be obtained analogously. 

Note that the vector (m + 2 — A, — e q , —e q , . . . , — e q ) 
of the first row of the matrix on the rhs of Eq. 
can be looked upon as the sum (i.e., a linear combina- 
tion with all scalars being 1) of the following to + 3 vec- 
tors: (m + 2 - A, q , q , q ), (0, -e g , q , . . . , 9 ), and 
(0, q , -e q , ...,0 q ), (0,0 q ,0 q ,..., -e q ), where 0, is 
a zero vector with order N g — 1. For all these to + 3 
vectors with identical order N g +i, only one entry is non- 
zero (to + 2 — A or— 1), while the N g+ i — 1 entries are 
zeros. According to the properties of determinants, we 
can rewrite P 3 +i(A) as a liner function of the first row 
on the rhs of Eq. ([551) when the remaining rows are held 
fixed. Concretely, P g+ i(A) can be expressed as 



A 0„ 



— eJ 



q 

o 



O Q 9 

o o 

-e„ • 



9 

o q 

ej O Q 9 



O 



o o 



o 
o 



Q fl 



u 9 

o 
o 





-e 1 
— e ' 



— e ' 



-e, q 

Qg O 

O Q s 



o o 



o q o q 



O Q„ 



O O 



"9 

o 
o 



Q 



g 

-e q 
O 
O 



(Al) 



It is easy to see that the last m + 2 determinants on 
the rhs of Eq. (|A1[) are equal to one another. We use 
Pj+i(A) to denote the first determinant and use Pj+\(A) 
to represent any of the last m + 2 determinants. Thus, 



the matrix corresponding to the determinant, we have 



P 3+1 (A)=P« 1 (A) + (to + 2)P^ 1 (A) 



(2) 



(A2) 



Below we will calculate these two quantities P^\(A) and 
^ ( +\(A). 

For Pg^\(A), it is obvious to have 

P 9 ( i\(A) = (m + 2-A)[Q 9 (A)r+ 2 . (A3) 
With regard to P^(A), after some transformations of 




[Q 9 (A)r +1 [|L 9 l + (A-l)|Q s |] 
[Q 9 (A)] m + 1 [P ff (A) + (A-l)0 g (A)].(A4) 
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Plugging Eqs. (K3\i and (1X4)) into Eq. (lA"2|) . we obtain 

Pg+l 



where r is a zero vector with order N g — 2. 

Appendix B: Using recurrence relations for relevant 
polynomials to determine EMFPT of Vicsek fractals 

In order to demonstrate the universality of our method 
for calculating EMFPT of self-similar tree-like fractals, 
here we apply our proposed method to compute the 
EMFPT for the Vicsek fractals jilli^. which is one of 
the most important and frequently studied regular fractal 
classes. 

The Vicsek fractals are built in an iterative way [HIHH 
controlled by the two parameters / and g. Denote by V/ jS 
(/ > 2, g > 1) the Vicsek fractals after g iterations (gen- 
erations). The initial construction (g = 1) is a starlike 
cluster composed of / + 1 nodes arranged in a crosswise 
pattern with / peripheral nodes connected to a central 
node. This corresponds to V/,i. Hereafter, we call those 



which is exactly Eq. (|43| in the main text. 



I 

nodes farthest from the central node at any iteration as 
outmost nodes. For g > 2, V/ )fl is obtained from V/ )S _i. 
To obtain V/ jS , we generate / identical copies of Vy j9 _i 

(denoted by yf g _ 1} i = 1,2,...,/) and arrange them 
around the periphery of the original V/ iS _i (denoted by 

y^g-i), then we add / new edges, each of them connects 
an outmost node in one of the / corner copy structures 
and that of the original central structure as shown in 
Fig. El 

After introducing the Vicsek fractals, in what fol- 
lows, we will apply our method to determine exactly the 
EMFPT between two nodes over the whole Vicsek fractal 
family. Note that in the case without confusion, we will 
use the same notations as that used above for the tree- 
like fractals T g . It should be mentioned that since the 
number of recurrence relations for related polynomials 
increases with parameter /, below we only provide the 
calculation detail of the EMFPT for a particular case of 



v) = (m + 2-\)[Q g (\)} 



m+2 



+(m + 2)[Q 9 (A)r +i [P g (A) + (A — l)Q g {\)] 



(to + 2)[Q g {\)] m+1 P g {\) + (to + l)A[Q g (A)] T CA%id R g +i{\) are as follows: 



In a similar way, we can easily derive the recursion 
relations for Q g +i(A) and i? g+ i(A) expressed in Eqs. (|44|) 
and P5)l . The concrete derivation processes for Q g+ i(X) 



(A) 



[Q g wr +1 



[QgW] 



m+l 



m + 2- X 
-e T 

e r 

TO + 1 r 

-el R, 



R„ 



(m + l)R g (\)[Q g (\)Y 







1 — A — e r 

+ (to + l)R g (X)[Q g (X)r[P g (X) + (A — l)Q g (X)} 
= [Q 9 WY n+1 [(m + l)|R fl | + |Q ff |] + (m + l)R g (X)[Q g (X)] m [P g (X) + (A - l)Q g (X)] 
- [Q 9 (A)P +1 [(m + l)R g (X) + Q g (X)] + (to + l)R g (X)[Q g (X)] m [P g (X) + (A - l)Q g (X)] 
= [Q a W] m+2 + (m + l)XR g (X)[Q g (X)} m+1 + (to + l)R g (X)[Q g (X)] m P g (X). 



(A6) 



Rg+lW = [QgM] 



[Q B W] 



+ 2 — A — e, r — e r 
ej R 3 O 
O R„ 

0, 

o 

R, 



1 - A 

-e, T 
-el 



Rg 

O 



+m[R g (X)] 2 [Q g (X)} 
[Q g (X)] m (R g (X) 



m—1 \ 



9 

1 - A -e r 
— e r R 9 
2r/n /wim-l 



m{R g (X)f{Q g (X)Y 




-e] 



Q 9 



-ej 



0, 

Rg 

o 



O 

R„ 



(A-l)Qg(A)] 



R g (X) 



1 - A 
-ej 



-e, 

R n 



Jig (A) 



TO + A 

-ej 



0, 
R„ 



+m[R g (X)] 2 [Q g (X)r- 1 [P g (X) + (A - l)Q g (X)} 
[Q g (X)] m (2R g (X)\Q g \ + (to + A)i?g(A)|Rg|) + TO[ J Rg(A)] 2 [Q s (A)] m - 1 [P 9 (A) 

TO[i?g(A)] 2 [Qg(A)]™- 1 [Pg(A) 

m[RJX)] 2 [QJX) 



[Q g (X)} m {2R g (X)Q g (X) + (to + X)[R g (X)} 2 } 



(A 
(A- 



1)Q S (A)] 
1)Q 9 (A)] 



2i?„(A)[gg(A)] m+i + (to+ l)X[RJX)nQJX)Y 



^(A), 



(A7) 



11 



FIG. 5: The first several iterative processes of a particular 
Vicsek fractal 1/4,3. The open circles express the starting 
structure V4 1 . 



/ = 3. 

For the convenience of description, we define V g (X) as 
the characteristic polynomial of the Laplacian matrix L 9 
corresponding to ¥3.3, i.e., 



V g (X) = dct(L s -AI 9 



(Bl) 



In order to obtain the sum of the reciprocal of all nonzero 
eigenvalues of L s , viz., all nonzero roots of polynomial 
V g (A) , we denote V g (X), V" 9 2 (A), and V g 3 (A) as deter- 



where V q {X) 



minants of sub- matrices of (L g — AI S ), 
(i = 1,2,3) is obtained by deleting from (L s — XI g ) i 
rows and columns corresponding to i outmost nodes that 
belong to i different copy structures forming V3. g , i.e., 



3,9- 



V 



(2) 

3,9-1' 



3,9- 



The polynomials V g (X), V g 1 (X), V^A), and F g 3 (A) obey 
the following recursive relations: 

V a+1 (X) = [V g (X)} i + [6V Q 1 (X)+3V"(X) + V Q 3 (X)}[V g (X)f 



+ [9[V g l (X)} 2 +3Vg l (X)V g 2 (X)}[V g (X)} 2 

+4lV g 1 (X)fV g W, 
Vj+iW = [Vg 1 (X) + V^X)][V g (X)f + [5{V g 1 (X)f 



(B2) 



and 



V g 3 +1 (X) = [V q 1 (X) + V q 2 (X)} 3 V g (X) + 3[V q \X)} 4 



-9[V g \X)]XW 

'[Vg\X)]XW- 



(B5) 



Using above recursion relations one can compute the 
sum of the reciprocal of the non-zero roots of V g (X). To 
this end, we define another polynomial V^'(A) as 



K(x) 



Vg(X). 



(B6) 



Obviously, the sum of roots of V^'(A) is equal to the sum of 
non-zero roots of V g (X). According to the above-obtained 
results, we have 



V,(o)' 



(B7) 



where v' g (0) is the constant term of V g (X) and v'(l) is 
the coefficient of term A of V^'(A). Then, the problem 
of determining (H) g is reduced to finding the quantities 
v' g (0) and v' g (l). 

From Eqs. (|B2[) - f|B5[) we can easily get the following 
four recursive equations: 

y 3 ' +1 (A) = X 3 [V g '(X)] i + X 2 [6V g 1 {X) + 3V g 2 {X) + V g a (X)] 

[K(A)] 3 + X[9[V g \X)} 2 + 3FVA)V 2 (A)][K'(A)] 



+4[^(A)] 3 ^(A), 
V^X) = A 3 [U 9 1 (A)+y 3 2 (A)][^(A)] 3 + A 2 [5[^ 1 (A)] 



(B8) 



v g 2 +1 W 



+7V g 1 (X)V g 2 (X) + [V g 2 (X)} 2 + ^(A)y ff 3 (A)][y;(A)] 2 

+Xl5[V g \X)f + 5[^(A)] 2 I/ 2 (A)]T/'(A) 

+ [^(A)] 4 , (B9) 

A 2 [^(A) + lf(A)] 2 [V;'(A)] 2 + A[4[^(A)] 3 
+9[V}{X)fV 2 {X)+W q 1 {X)[V q 2 {X)f 



-[V g 1 {X)fV^{X)]V g , {X) 
-3[^(A)] 3 V; 2 (A), 



^(A)] 



(BIO) 



and 



V g 3 +1 (X) 



A[K 1 (A) + U 2 (A)] 3 K'(A)+3[K 1 (A)r 



+7V g 1 (X)V g 2 (X) + [Vf(A)] 2 + ^(A)V S 3 (A)][V S (A)] 2 
+[5{V g 1 (X)} 3 + 5[V g 1 (X)] 2 V g 2 (X)]V g (X) 

+[n 1 (A)] 4 , 



-9[V g \X)]XW 
-[V g \X)fV*(X). 



9 wr 

(Bll) 



^;+i(a) 



= [VgHX) + V 2 (X)] 2 [V g (X)} 2 + [4[Vg\X)] s 
+9[K 1 (A)] 2 V 2 (A) + 4V 1 (A)[K 2 (A)] 2 



+ [V g 1 (X)] 2 V g 3 (X)]V g (X) 
+3[^ 1 (A)] 3 V g 2 (A), 



(B3$ased on these relations we can find the values for v q (0) 
and v' g {l). 

Let us first evaluate v' g (0). To this end, we use v g (0), 
Vg(0), and v g (Q) to denote separately the constant terms 
of V}(\), V 2 (X) and V*(\). Making use of Eqs. (fB8|) - 
(jBlip . we can derive the following important relations: 



(B4) 



<+i(0)=4[^(0)]V(0). 



(B12) 
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^+i(0) = K(0)] 4 , (B13) 



« 2 +1 (0) = 2[^(0)] 4 + 3[«*(0)] 3 w*(0) , (B14) 

and 

^ 3 +i(0) = 3[^(0)] 4 + 9[^(0)] 3 ^(0) 

+6[^(0)] 2 b g 2 (0)] 2 + [^(0)]\ 3 (0)(B15) 

Considering ^( 0) = -4, ^(0) = 1, w 2 (0) = 2, and 
vf(0) = 3, Eqs. (|B12[) - (|B15I) are resolved to obtain 

<(0) = -4* , (B16) 



w£(0) = l, (B17) 



V 2 (0) = 33-l, (B18) 

and 

u 3 (0) = ^(3*-l) 2 . (B19) 

With these obtained results, we continue to calculate 
v'(l). To achieve this goal, we define Vg(l) as the co- 
efficient of term A of polynomial V g (X). According to 
Eqs. (|B8l) - (|Blip . the quantities v' g (l) and v g (l) can be 
represented recursively as follows: 

= 4[^(0)]V fl (l) + 12[^(0)]X(0)^(l) 
+9[^(0)] 2 K(0)] 2 + 3^(0)^0)^20) 
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